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Abstract. Based on self- consistent field theory we study a soliton gen- 
O ' eration in cw solid-state lasers with semiconductor saturable absorber. Var- 

O . ious soliton destabilizations, i.e. the switch from femtosecond to picosecond 

^ . generation ("picosecond collapse"), an automodulation regime, breakdown of 

soliton generation and hysteresis behavior, are predicted. It is shown that the 



>^| third-order dispersion reduces the region of the soliton existence and causes 



the pulse oscillation and strong frequency shift. 
OCIS codes: 140.7090, 140.3430, 190.5530 



1 Introduction 

Recent advances in ultrafast solid-state lasers have resulted in sub 6-fs pulses 
generated directly from the cavity of Ti: sapphire laser [1], which is close to 
the theoretical limit for optical pulses. The key role to mode-locking in this 
case plays the fast optical nonlinearities, such as electronic nonlinearities, 
responsible, in particular, for self-phase modulation (SPM) and self-focusing 

[2]. 

As was shown, the generation of extremely short pulses in the lasers with 
self-focusing or semiconductor saturable absorber [3] is possible due to forma- 
tion of a Schrodinger soliton. The relevant mechanism is the balance between 
SPM and group velocity dispersion (GVD), which reduces the pulse dura- 
tion and stabilizes it. This mechanism termed the soliton mode-locking [4] 
works together with Kerr-lensing or saturable absorption of semiconductor 
thus supporting the ultrashort pulse and preventing from the noise gener- 
ation. However, the presence of dissipative laser factors, such as saturable 
gain and loss, frequency filtering etc., complicates the soliton dynamics essen- 
tially, so that various nonstationary regimes in Kerr-lens mode-locked lasers 
[5, 6] and Q-switching instability in the cw lasers with semiconductor sat- 
urable absorber [7] are possible. Up to now a theoretical understanding of 
the above-mentioned issues is lacking and challenges for investigators efforts, 
that can be of particular interest for the for generation dynamics and pulse 
characteristics control in femtosecond region. 

Here we present the results of our study of the soliton generation in cw 
solid-state lasers with semiconductor saturable absorber. The main focus 
of investigation was set on the peculiarities of the transition between fem- 
tosecond (soliton) and picosecond generation. We shall demonstrate that 
such a transition is accompanied by the threshold and hysteresis phenom- 
ena. Based on soliton perturbation theory, the numerical simulations of two 
different experimental situations have been performed. The first one corre- 
sponds to the variation of control parameters (dispersion or pump power), 
when for every new value of control parameter the laser is turned on afresh. 
The second situation is for continuous variation of control parameter dur- 
ing a single generation session. We demonstrate also that the third order 
dispersion destabilizes the soliton and produces a strong frequency shift of 
generation wavelength with respect to the gain center. 



2 Model 



We analyzed the field evolution in the distributed laser system containing sat- 
urable quasi-four level gain crystal, frequency filter, GVD and SPM elements, 
two level saturable absorber and linear loss. We assumed, that in noncoherent 
approximation the saturable absorber action can be described by operator 

~'^'^uWdt 1^1' ■^here 7 is the initial saturable loss, E = J \a{tf)\ dti is 

—00 
the energy fiuency passed through the absorber to moment t, t is the local 

time, Es is the loss saturation energy fiuency; the expansion of the denom- 
inator into the time series in ^^ is further supposed. The gain saturation is 
due to the full pulse energy. The ratio of the loss saturation energy to the 
gain saturation energies r under assumption of the equal cross sections of 
the generation mode at semiconductor modulator and Ti: sapphire active 
medium is 1.25 x 10"'^ for Es = 100 jiJ /cm? (compare this with the figure 
for semiconductor saturation energy reported in [3]). 

The evolution of generation field a obeys the following operator equation: 
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were z is the longitude coordinate normalized to the cavity length, or 
transit number, a is the saturated gain, / is the linear loss, A;2,3 are the 
second- and third-order dispersion coefficients, respectively, the energy fiu- 
ency is normalized to Eg . The time is normalized to the inverse bandwidth of 
the absorber 1/ta . With this normalization the SPM coefficient f3 is ^^^^f^, 
where L is the length of crystal, n and n2 are the linear and nonlinear coeffi- 
cients of refractivity, respectively, A is the generation wavelength. The term 
in square parenthesis stands for the frequency filtering. The transmission 
band of the filter and the absorption band of semiconductor were assumed 
to coincide. An expansion of Eq. (1) up to second-order in d/dt results in 
nonlinear Landau-Ginzburg equation [4, 8], which we do not write out here 
due to its complexity. 

Since an exact general solution of Eq. (1) is unknown, we sought for the 
approximated quasi-soliton solution in the form a(z,t) = Qo '^^'tfu ^w-!"\i+'^ ■• 
where a^ is the amplitude, tp is the width, oj is the frequency mismatch from 
filter band center, \i/ is the chirp, ip and 5 are the phase and time delays 
after the full cavity round trip, respectively. In the frame of aberrationless 



approach [9], the substitution of this solution in Eq. (1) with following 
expansion in t yields the following set of ordinary differential equations for 
evolution of pulse parameters: 
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where v = 1/tp is the inverse pulse width, and the solutions for time and 
phase delays: 

5 = [al (cj2 + ^;2 _ 1) _ 2v^iu^)] 7 exp (-^) - 2iu'^v^ - 2k2iuv^+ 
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Assuming that pulse width is much shorter than the cavity period Tcav , 
the equation for the gain evolution is as follows: 

da , . ^ 2aTal a , , 

— = (a„-a)P-— -^--, (3) 

where P = ^"J '' is the dimensionless pump intensity, a is the absorption 
cross-section at the pump wavelength, hv is the energy of pump photon, a^ 
is the maximal gain at the full population inversion. P = 10~^ corresponds 
to the pump power of IW for 100 /xm diameter of pump mode. 



The steady-state solution of Eqs. (2, 3) describes the stable ultrashort 
pulses generation. The system (2, 3) was solved by the forward Euler method 
with number of iterations 10^ and the accuracy of 10"^ . 

3 Discussion 

First, we shall study the situation when for each new value of control param- 
eter the generation is formed starting from the noise spike as initial approx- 
imation. 

The normalized width of the stable pulse versus GVD coefficient is pre- 
sented in Fig. 1 for different amount of SPM. It is well known that in 
solid-state lasers contrary to the dye lasers a soliton mode locking with a 
slow saturable absorber is impossible because of small contribution of dy- 
namic gain saturation (a small r in our notations) but at the same time the 
presence of SPM can provide the stable soliton mode locking [4]. As was 
expected, our calculations confirmed this conclusion for the situation with 
no SPM present in the system. However, with some SPM another pulse 
stabilizing mechanism comes to play a role. Introducing the SPM stabilizes 
the picosecond generation (region A in Fig. 1) due to the action of negative 
feed-back [8] (top-left segment of curve 1 in Fig. 2) . The mechanism of this 
feed-back is as follows: an increased pulse intensity causes a stronger chirp- 
ing, i.e. a wider pulse spectrum and, consequently, a higher loss at bandwidth 
limiting element and, reverse, a decreased intensity minimizes the chirp and 
the spectral width, reducing the loss at spectral filter. 

There is a section at the curves 1-4 (Fig. 1, 2) where the Schrodinger 
soliton exists within limited window of negative dispersion. The pulse du- 
rations lie in femtosecond region with the minimum close to zero GVD, the 
pulse intensities are much higher and the chirp is very small (central part 
of curve 1 in Fig. 2). Calculation showed that contrary to ps-case and the 
situation with non-zero third-order dispersion the soliton has no frequency 
shift from filter band center and its energy is equal to critical energy for 
first-order Schrodinger soliton E = -j-|i . 

At some negative GVD the switch from femtosecond to picosecond regime 
("picosecond collapse" marked by arrow, see curve 1 in Fig. 1) takes place. 
Formally, these two types of ultrashort generation correspond to domination 
of either nondissipative (case of Schrodinger soliton generation) or dissipative 
(case of picosecond generation) terms in Eq. 1, respectively. 

For the stronger SPM (curve 2 in Fig. 1) the interval of GVD, where 
the femtosecond generation takes place, broadens which is accompanied by 
the shortening of the pulse width (compare curves 1 and 2 in Figs. 1,2). 



As is seen, the ps-collapse occurs now at much higher (both in negative and 
positive region) dispersions. 

Further increasing of SPM transforms the character of sohton destabiUza- 
tion (curve 3 in Fig. 1). In this case, the femtosecond generation is possible 
only near zero dispersion; the higher negative GVD destroys the soliton-like 
pulse (discontinuity of curve 3). Left to the region where no quasi-soliton 
solution for Eq. 1 exists, there appears the region of the pulse automodu- 
lational instability (region B). Here, the pulse parameters oscillate with the 
period of ~ 100 cavity round-trips (~ 1 /is ) as is illustrated in Fig. 3. The 
oscillating pulse has quasi-Schrodinger soliton type with a small chirp (region 
B in Fig. 2) and frequency shift. 

Further increase of SPM (curve 4 in Fig. 1) reduces the interval of quasi- 
soliton existence. Comparing curves 1 - 4 of Fig. 1 one may conclude that 
the optimal amount of SPM exists providing the femtosecond generation in 
the widest interval of GVD. 

Let's consider the pump intensity P - another control parameter of the 
system switching the generation between ps- and fs-regimes. As is seen from 
Fig. 4 (triangles and arrow 1), the increase of the pump switches the laser 
from ps- to fs-generation. This occurs when the pump intensity and, con- 
sequently, the generation pulse intensity becomes high enough for SPM to 
compensate the GVD. 

We can investigate the system behavior depending on control parameters 
(pump, GVD) by two different cases. First one is when we measure the pulse 
characteristics at the some value of control parameters then turn off the laser 
where upon turn on the laser and iterate this process with another value of 
control parameter. Mathematically it corresponds to receiving of the pulse 
as stationary solution of Eq. 1 from the noise spike. This case is depicted at 
Figs. 1, 2 and triangles in Fig. 4. Second case is when one receives ultrashort 
pulse at the definite value of the control parameter and then begins smooth 
variation of this parameter. This case is depicted at Fig. 4, b and Fig. 5. 

As is seen from Fig. 4, the dependence of the pulse duration on the pump 
intensity demonstrates a hysteresis behavior. The pump variation from the 
small values to the large and backward produces the following sequence of 
generation regimes (trajectory 2 in Fig. 4): the pump increase from 10~^ to 
2 X 10""^ does not bring the system from ps- to fs-generation, however the 
subsequent pump decrease from 2 x 10"*^ to 1.67 x 10~^ switches the system 
to fs-regime. The range of fs-generation in terms of P (solid curve) is wider 
than for the case denoted by triangles, which was plotted analogous to curves 
of Figs. 1 and 2. If the trajectory starts from large P, the picosecond regime 
does not set and the femtosecond soliton forms at P ~ 1.67 x 10""' (solid 
curve). In this case, the interval of the soliton existence is wider than for 
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regime denoted by triangles. 

The system demonstrates the hysteresis behavior under the variation of 
GVD, too (Fig. 5): changing the GVD from the negative values to zero 
and than to the positive values does not cause the switch from ps- to fs- 
generation. On the other hand when moving from the positive values of 
dispersion to negative and backwards one can see the dramatic change in 
system's behavior. The decrease of GVD starting from the positive values, 
where the ps-generation takes place, causes an abrupt switch to fs-generation 
near zero GVD. Comparing curve 2 in Fig. 1 and the curve in Fig. 5 that 
were plotted for the same parameters but for different physical situations 
(described above), one may conclude that the region of the soliton existence 
is wider in term of GVD for the latter case (when we alter the GVD) than 
for the former (when generation is formed from the noise spike). The switch 
from fs to ps generation under inverse movement along the GVD axis takes 
place at another value of GVD thus producing a hysteresis feature (see Fig. 
5). 

Since the Schrodinger soliton formation occurs at small negative GVD, 
the contribution of third-order dispersion k^ in this region may be essential. 
The presence oi k^ ^ gives rise to the following effects: 1) the strong pulse 
frequency shift which depends on the sign of k^ (compare curves 2, 3 and 
4 and curve 1 in Fig. 6) arises; 2) the region of fs-generation narrows; 3) 
an additional destabilizing factor, the frequency shift oscillations, appears 
(regions B for curves 2, 3, 4). 

4 Conclusion 

In summary, based on the self-consistent field theory we have investigated 
the characteristics of Schrodinger soliton in cw solid-state laser with semicon- 
ductor saturable absorber. We demonstrated numerically, that the formation 
of soliton has threshold character, so together with first (free running) and 
second (mode-locking) threshold the threshold of femtosecond generation ex- 
ists. Three main destabilization scenarios were demonstrated: the switch to 
picosecond generation (picosecond collapse), the switch to automodulation 
mode and the breakdown of soliton-like pulse. The switch between pico- 
and femtosecond generation has the hysteresis character. The contribution 
of high-order GVD reduces the interval of quasi-soliton existence, gives rise 
to the frequency shift and the pulse destabilization. 

Our results may be useful for design of self-starting cw solid-state lasers 
with controllable ultrashort pulse characteristics. 

This work was supported by National Foundation for basic researches 



(grant F97-256). 
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Figure 1: Pulse width tp versus GVD coefficient ^2 (time is normalized to 
ta ). Region A, picosecond generation, region B, autooscillation regime for 



parameters of curve 3. P 

Tcav = 10 ns, I = 0.01, 7 = 
0.1 (4). 



= 5 X 10 , am = 0.5 , ta = 2.5 fs, T^ = 3 yus , 
0.05 ,k3 = 0.p = 0.001 (1), 0.01 (2), 0.05 (3), 
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Figure 2: Chirp \1/ versus GVD coefficient k2 . The parameters correspond 
to Fig. 1. Curve 3 and 4 exist under k2 > —3,-1, respectively. 
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Figure 

k2 = -20 



3: Peak intensity / normalized to Est^^ ^ versus round-trip number z. 
^n ^g = , P = 5 X 10"^ , p = 0.05. 
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Figure 4: Pulse width tp versus pump intensity P. Arrow 1 - switch be- 



tween ps- and fs-generation (triangles), 
parameters correspond to Fig. 1. 
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Figure 5: Pulse width tp versus GVD coefficient /c2 . Parameters correspond 
to curve 2 in Fig.l. 
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Figure 6: Frequency shift lu versus GVD coefficient k2 ■ Parameters corre- 
spond to curve 3 in Fig. 1. h = (1), -5 (2), -10 (3), 10 (4). The initial 
conditions correspond to Figs. 1, 2. 
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